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1 Introduction 



For a long time it was widely believed that the admixture of the pairs of 
strange quarks in the nucleons is small. The main justification of this picture 
was the constituent quark model where there is no room for strange quark in 
the nucleon. It has been known for a while that this picture is not quite true: 
In scalar and pseudoscalar channels one can expect a noticeable deviation 
from this naive prediction. This is because, these channels are very unique 
in a sense that they are tightly connected to the QCD-vacuum fluctuations 
with 0"*", 0~ singlet quantum numbers. Manifestation of the uniqueness can 
be seen, in particular, in the existence of the axial anomaly (0~ channel) 
and the trace anomaly (0"'" channel). Nontrivial QCD vacuum structure tells 
us that one could expect some unusual properties when we deal with those 
quantum numbers. 

As we now know , this is indeed the case. In particular, we know that the 
strange quark matrix element (A^|ss|A^) does not vanish and has the same 
order of magnitude as {N\dd\N). This information can be obtained from the 
analysis of the so-called a term Similarly the analysis of the " proton 

spin crisis" essentially teaches us that the spin which is carried by the strange 
quark in the nucleon is not small as naively one could expect, see e.g. the 
recent review 0. 

Another phenomenological manifestation of the same kind is the very old 
observation that in the scalar and pseudoscalar channels the Zweig rule is 
badly broken and there is substantial admixture of s quarks in the scalar 
mesons /o(980) (was S*), ao(980) (was 6), /o(1300) (was e), as well as 
in the pseudoscalar mesons r] and rj'. At the same time, in the vector 
channel the Zweig rule works well. Phenomenologically it is evident in 
e.g. the smallness of the (p — u mixing. In terms of QCD such a small- 
ness corresponds to the numerical suppression of the nondiagonal correlation 
function / dx{0\T{s'y^s{x),U'y^u{0)}\0) in comparison with the diagonal one 
/ dx{0\T{u'-y^u{x),U'-y^u{0)}\0) . In the scalar and pseudoscalar channels di- 
agonal and non-diagonal channels have the same order of magnitude. We 
believe that analysis of such kind of the correlation functions is an appropri- 
ate method for a QCD- based explanation of the unusual hadronic properties 
mentioned above. 

In this talk we present some general methods and ideas for the analysis 
of the nucleon matrix elements from a non-valence operator. The ideology 
and methods (unitarity, dispersion relations, duality, low-energy theorems) 
we use are motivated by QCD sum rules. However we do not use the QCD 
sum rules in the common sense. Instead, we reduce one complicated prob- 
lem (the calculation of non-valence nucleon matrix elements) to another one 
(the behavior of some vacuum correlation functions at low momentum trans- 
fer). One could think that such a reducing of one problem to another one 
(may be even more complicated) does not improve our understanding of the 
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phenomenon. However, this is not quite true: The analysis of the vacuum 
correlation functions with vacuum quantum numbers, certainly, is a very dif- 
ficult problem. However some nonperturbative information based on the low 
energy theorems is available for such a correlation function. This gives some 
chance to estimate some interesting quantities. 

2 Strangeness in the nucleon, 0+ channel. 
2.1 First estimations 

We start by calculating the strange scalar matrix elements over the nucleon, 
assuming an octet nature of SU{3) symmetry breaking. We follow to ref.[Q( 
see also the book |^ for a review) in our calculations 0, but with a small 
difference in details. We present these results for completeness of the talk. 

The results of the fit to the data on nN scattering presented in [@] lead 
to the following estimates for the so-called a term 

!!hL±I!^ (^p\uu + dd\p) = (64 ± 8MeV) . (1) 

(Here and in what follows we omit kinematical structure like pp in expressions 
for matrix elements.) Taking the values of quark masses to be rriu = 5.1 ± 
Q.9MeV , TUd = 9.3 ± lAMeV , m, = 175 ± 25MeV 0, from (D we have 

{p\uu + dd\p) ~ 9. (2) 

Further, assuming octet-type SU{?)) breaking to be responsible for the mass 
splitting in the baryon octet, we find 

(j)\uu - dd\p) = "^^'"^^ = 0.7, (3) 

{p\uu + dd- 2ss\p) = ^ 3_4_ a) 

Here ms, "^2; are masses of H, S, A hyperons respectively. The values (^, 
(^ are quite reasonable: the former is close to the difference of the number 
of u and d quarks in a proton (should be 1), and the latter is close to the 
total number of valence quarks u and dm. a. nucleon (see below), from (@-|^) 
one obtains: 

{p\uu\p) ~ 4.8, (5) 
{p\dd\p) ~ 4.1, (6) 
{p\ss\p) - 2-8. (7) 
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We should mention that the accuracy of these equations is not very high. 
For example, the error in the value of the a term already leads to an error 
of order of one in each matrix element discussed above. However, these very 
simple calculations explicitly demonstrate that the strange matrix element 
is by no means small. 

We would like to rewrite the relations (0-0) to separate the vacuum con- 
tribution to the nucleon matrix element from the valence contribution. In 
order to do so, let us define 

{p\qq\p) = {p\qq\p)o + {p\qq\p)i, (8) 

where index labels a (sea) vacuum contribution and index 1 a valence 
contribution for a quark q. We assume that the vacuum contribution which 
is related to the sea quarks is the same for all light quarks u, d, s. Thus, the 
nonzero magnitude for the strange matrix elements comes exclusively from 
the vacuum fluctuations. At the same time, the matrix elements related to 
the valence contributions are equal to 

~ (4.8 - 2.8) ~ 2, (9) 

{p\dd\p)i ~ (4.1 - 2.8) ~ 1.3. (10) 

These values are in remarkable agreement with the numbers 2 and 1, which 
one could expect from the naive picture of non-relativistic constituent quark 
model. In spite of the very rough estimations presented above, we believe we 
convinced a reader that : 

a) a magnitude of the nucleon matrix element for ss is not small; 

b) the large magnitude for this matrix element is due to the nontrivial QCD 
vacuum structure where vacuum expectation values of u, d, s quarks are de- 
veloped and they are almost the same in magnitude: (0|(i(i|0) ~ (OImmIO) ~ 
{0\-ss\0). 

Once we realized that the phenomenon under discussion is related to 
the nontrivial vacuum structure, it is clear that the best way to understand 
such a phenomenon is to use some method where QCD vacuum fluctuations 
and hadronic properties are strongly interrelated. We believe, that the most 
powerful analytical nonperturbative method which exhibits these features is 
the QCD sum rules approach 

In what follows we use the QCD sum rule method in order to relate 
hadronic matrix elements and vacuum characteristics. Let me emphasize 
from the very beginning that we do not use the QCD sum rules in the stan- 
dard way: we do not fit them to extract any information about lowest reso- 
nance (as people usually do in this approach), we do not use any numerical 
approximation or implicit assumptions about higher states. Instead, we con- 
centrate on the qualitative relations between hadronic properties and QCD 
vacuum structure. We try to explain in qualitative way some magnitudes 
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for the nucleoli matrix elements which may look very unexpected from the 
naive point of view. At the same time those matrix elements can be easily 
understood in terms of the QCD vacuum structure. 

We close this section with the formulation of the following question: 
Q: What is the QCD explanation of the unusual properties mentioned above? 
(in particular, the large magnitude for the strange nucleon matrix element, 
a special role of the scalar and pseudoscalar channels et cetera). Our answer 
on this question is: 

A: Hadronic matrix elements with 0^ quantum numbers are singled out be- 
cause of the special role they play in the QCD vacuum structure. The next 
section switches this answer from a qualitative remark into the quantitative 
description. 

2.2 Strangeness in the nucleon and vacuum structure 

To study the problem of calculation (A^|ss|A^) using the QCD -sum rules 
approach, we consider the following vacuum correlation function 0]: 

T(g2) = j e''^^dxdy{Q\T{r^{x)My).mm (H) 

at — oo. Here rj is an arbitrary current with nucleon quantum num- 
bers. In particular, this current may be chosen in the standard form rj = 
e°'^'^'-)^d'^{v!'C'^^u'^). Note, however, that the results obtained below do not 
imply such a concretization. For the future convenience we consider the unit 
matrix kinematical structure in ([TT|). 

This is the standard first step of any calculation of such a kind: Instead 
of direct calculation of a matrix element, we reduce the problem to the com- 
putation of some correlation function. As the next step, we use the duality 
and dispersion relations to relate a physical matrix element to the QCD- 
based formula for the corresponding correlation function. This is essentially 
the basic idea of the QCD sum rules. 

In our specific case ( pH]) due to the absence of the s -quark field in 
the nucleon current ?], any substantial contribution to T[q^) is connected 
only with non-perturbative, so-called induced vacuum condensates, see Fig.l. 
Such a contribution arises from the region, when some distances are large: 
iy — 0)^ ~ (y — x)"^ ^ (a; — 0)^. Thus, it can not be directly calculated in 
perturbative theory, instead we code the corresponding large-distance infor- 
mation in the form of a bilocal operator 

K = ij dy{Q\T{ss{y),uu{Q)m. (12) 

see Fig.l The similar contributions were considered at first time in ref. ||1CI||. 
For the different applications of this approach when bilocal operators play 
essential role, see also refs.[|n|]. 
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Along with consideration of the three-point correlation function ([TlD, we 
would like to consider the standard two-point correlator 

P(g2) = I e''^^dx{0\T{vix),r](0)}\0), (13) 

see Fig. 2. The correlator (|13]) is determined by the nucleon residues {0\ri\N) 
and some duality interval 5*0. At the same time the correlator ([TT|) includes 
the information on the nucleon matrix element (A^|ss|A^) also. Comparing 
(pj]) with (|13D at — — »• cx), we arrive to the following relation ||^: 

— in 

{N\ss\N) ~ J-^K, (14) 

where m is the nucleon mass. The main assumptions which have been made 
in the derivation of this relation are the following. First, we made the stan- 
dard assumption about local duality for the nucleon. In different words we 
assumed that a nucleon saturates both correlation functions with duality in- 
terval 5*0. The second assumption is that the typical scales (or what is the 
same, duality intervals in the limit —q^ — > oo in the corresponding sum rules 
( (|lT]) and (|l^) ) are not much different in magnitude from each other. In 
this case the dependence on residues (0|?7|iV) is canceled out in the ration 
and we are left with the matrix element (A^|ss|A^) (|l^ we are interested in. 

Note, that both these assumptions are very likely to be satisfied because 
we know that in most cases the lowest state (nucleon) does saturate the sum 
rules. If it does, than the typical scale (which in variety of sum rules is one 
and the same and of order of IGeV"^) guarantees that the duality intervals are 
likely to be very close to each other. Anyway, the quantitative analysis of the 
corresponding sum rules is possible, however it is not our main goal; rather 
we want to demonstrate the relation between matrix elements like (A^|ss|A^) 
and the corresponding vacuum properties which are hidden in the correlator 
K (|12D. In principle one could analyze the sensitivity of the corresponding 
QCD sum rules to the lowest state, nucleon. Once it is demonstrated, we 
believe that the accuracy of our formula (0) is of order 20% — 30% which is 
a typical error for the sum rule approach. 

Thus, the calculating of (A^|ss|A^) reduces to the evaluation of the vac- 
uum correlator K. Fortunately, sufficient information about the latter comes 
from the low -energy theorems. We note also that this method of reduc- 
ing the nucleon matrix elements to that of the vacuum correlator is directly 
generalized to cover the arbitrary scalar O5 or pseudoscalar Op operator^: 

—mNN r 

{N\Os\N) ^-^^ J dy{0\T{Os,uuiO)m, (15) 

^We assume of course that these operators do not contain u,d quarks. Otherwise an 
additional contribution which conies from the smaU distances must be also included. 
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{N\Op\N) ^ "^/^'^^^ z / dy{0\T{Op,mj,u{0)m- (16) 



The estimation of the nonperturbative correlator K can be done by using 
some low-energy theorems. In this case K is expressed in terms of some 
vacuum condensates 0: 

K = tJ dy{0\T{ss{y),uuiOm ^ yvi^ ^ OMGeV, (17) 

\ 71- flu/ 

where b = yA^c ~ |^/ = 9 and we use the standard values for the vacuum 
condensates 0: 

(— ) ^ 1.2 ■ {qq) ~ -(250MeF)l 

71 

With the estimation (|17|) for K, our formula (|l^ gives the following expres- 
sion for the nucleon expectation value for ss 

, , , > 18 (qq) , , 

{v\-ss\v) ^ -m ■ YJ^^ ^ 2.4, (18) 

which is very close to the naive estimation (|^). Let us stress: we are not 
pretending to have made a reliable calculation of the matrix element {p\ss\p) 
here. Rather, we wanted to emphasize on the qualitative picture which 
demonstrates the close relation between nonvalence matrix elements and 
QCD vacuum structure. 

We close this section by noting that the method presented above gives 
very simple physical explanation of why the Zweig rule in the scalar and 
pseudoscalar channels is badly broken and at the same time, in the vec- 
tor channel the Zweig rule works well. In particular, the matrix element 
(A^|s7^s|iV) is expected to be very small as well as the corresponding cou- 
pling constant g^j^NN does. In terms of QCD such a smallness corresponds 
to the numerical suppression (10~^ — 10"'^) of the nondiagonal correlation 
function / dx{0\T{s'y^s{x) , U'y^u{0)}\0) in comparison with the diagonal one 
/ dx{0\T{u'y^u{x),U'-y^u{0)}\0) , see QCD-estimation in 0. In the scalar and 
pseudoscalar channels the diagonal and non-diagonal channels have the same 
order of magnitude. 

In the next few sections we discuss some applications of the obtained 
results. 



2.3 In the world where s quark is massless. 

We would like to look at formula ([T7|) from the different side. Namely, we 
note that K not only enters expression (0), but also determines the variation 
of the condensate (uu) with s quark mass: 

d r 

- — (uu) = -i / dy{0\T{ss{y),uu}\0) = -K ^ -QMGeV'^. (19) 
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To understand how large this number is and in order to make some rough 
estimations, we assume that this behavior can be extrapolated from physical 
value trLs ~ 175MeV till = 0. In this case we estimate that 



{uu) 



ms=175 

Such a decrease of | {uu) \ by a factor of two as varies from ~ llbMeV 
to = is a very important consequence of the previous discussions: Once 
we accept the relatively large magnitude for the nucleon matrix element 
{p\ss\p) ~ 2.4, we are forced to accept the relatively large variation of the 
light quark condensate as well. This statement is the direct consequence of 
QCD, see (|2D|). 

We note that this result does not seem very surprising since other vacuum 
condensates, e.g. {^C^^J) possess analogous properties [0. From the micro- 
scopic point of view, decrease of absolute values of vacuum matrix elements 
with the decrease of the s quark mass is expected since any topologically 
nontrivial vacuum configurations, e.g. instantons, are suppressed by light 
quarks. The corresponding numerical calculation is very difficult to perform, 
however a qualitative picture of the QCD vacuum structure definitely sup- 
ports this idea ||T3|. 



2.4 s quark and the nucleon mass. 

We would like to discuss here one more fundamental characteristic of the 
hadron world: the nucleon mass and its dependence on the strange quark. 
We start our discussion from the following well known result: the nucleon 
mass is determined by the trace of the energy -momentum tensor and 
in the chiral limit iJiu = rrid = mg = the nonzero result comes exclusively 
from the strong interacting gluon fields: 

h (y. 

^ = -^{N\—GUN), m„ = mrf = m, = 0. (21) 

However, as we know, in our world the strange quark is not massless, but 
rather it requires some (large enough) mass (~ 175MeV). As we have seen 
(pop, the nonzero mass of s quark considerably changes the vacuum properties 
of the world. Thus, we would expect that it might have strong influence on 
the nucleon mass as well. The main argument which supports this point of 
view is the same as before, and is based on our general philosophy that the 
nucleon matrix elements and vacuum properties are tightly related. So, if the 
strange quark has strong influence on the vacuum properties than its impact 
on the nucleon mass should also be strong. 

In order to check these reasons it would be useful to calculate the strange 
quark contribution into the nucleon mass directly and independently from 
the gluon contribution (pi]). Fortunately, it can be easily done by using our 
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previous estimation ([TSD for the nucleon matrix element and exact expres- 
sion for the trace of the energy-momentum tenzor with taking into account 
nonzero quark masses: 



b ,_,as 



m = +{N\Y.^MN) - ^im^GlJN), (22) 

q O 



where sum over q is sum over all light qurks u, d, s. One can easily see from 
(|l]) that u, d contribution into the nucleon mass does not exceed 7%; thus we 
can safely neglect this. At the same time, adopting the values ([7D, (p!8D for 



{p\ss\p) and ~ 175MeV 0, one can conclude that considerable part of 
the nucleon mass (about 45%) is due to the strange quark. In this case the 
gluon contribution into the nucleon mass is far away from the chiral SU{3) 
prediction (pTl) and approximately equals to 
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- ^{N\^GlJN) ~ 520MeV. (23) 



This rough estimation confirms our argumentation that a variation of the 
strange quark mass from its physical value to zero, may considerably change 
some vacuum characteristics as well as nucleon matrix elements. 

The simple consequence of this result is the observation that the quenched 
approximation in the lattice calculations is not justified simply because 
such a calculation clearly not accounting the fluctuations of the strange (non- 
valence) quark as well as vacuum fluctuations of u and d quarks. As we argued 
above, the nucleon mass undergoes some influence from s quark. 

How one can understand these results within the framework of the QCD 
sum rules? Let us recall that in the QCD sum rules approach an information 
about any dimensional parameter is contained in the vacuum condensates 
(uu), (G^^), .... As we discussed previously all these condensates varying with 
rris considerably. It is important that this variation certainly proceeds in the 
right direction: Absolute values of condensates decrease with decreasing m^. 
This leads to a smaller scale in the sum rules and finally, to the decrease of 
all dimensional parameters such as m. However, it is difficult to make any 
reliable calculations because of a large number of factors playing an essential 
role in such a calculation. 



2.5 Momentum distribution of the strangeness in the 
nucleon . 

We continue our study on the role of the strange quark in nucleon with the 
following remark. We found out earlier that the matrix element (A^|ss|A^) is 
not small; we interpreted this result as a result of strong vacuum fluctuations 
which penetrate into the nucleon matrix element. Now, we would like to ask 
the following question: What is the mean value of the momentum ( denoted 



8 



as {kY)s ) of the s quark inside of a nucleon? Let us note that this question 
is not a pure academic one. Rather, the answer on the question might be 
important for the construction of a more sophisticated quark model which 
would incorporate the strange context into the nucleon wave function. 

First of all, let us try to formulate this question in terms of QCD. We 
define the mean value (A;5.)s of the momentum carried by the strange quark 
in a nucleon by the following matrix element: 

{kl)mss\N) = {\N-s{i D^fs\N), (24) 

where i D^= i +gA'^^^ is the covariant derivative and is gluon field. 
The arrow shows the quark whose momentum is under discussion. 

We assume that the nucleon to be moving rapidly in the 2;-direction. 
We are interested in the momentum distribution in the direction which is 
perpendicular to its motion. Precisely this characteristic has a dynamical 
origin. Indeed, as we shall see in a moment, while we are studying a nucleon 
matrix element {k\)s, we are actually probing the QCD vacuum properties! 
The nucleon motion as a whole system with arbitrary velocity does not affect 
this characteristic. Thus, essentially, what we discuss is the, so-called, light 
cone wave function. Apart of the reasons mentioned above there are few 
more motives to do so: First of all, the light cone wave function {wf) with a 
minimal number of constituents is a good starting point. As is known such 
a function gives the parametrically leading contributions to hard exclusive 
processes. Higher Fock states are also well defined in this approach and can 
be considered separately. The second reason to work with a light cone wave 
function is there existence of the nice relation between that wf and structure 
function measured in the deep-inelastic scattering. We refer to the review 



paper [14] for the introduction into the subject. The relation to the standard 
quark model wave functions (see e.g.[|T3|) is also worked out. The relevant 
discussions can be found in ref. |jl6|. Besides of these, we have one more 
reason to work with the light cone wf: we believe that this is the direction 
where a valence quark model can be understood and formulated in the QCD 
-terms |]T7|. 

Anyhow, the formula ( p4D with the derivatives taken in the direction 



perpendicular to the nucleon momentum = {E,0±,pz), is very natural 
definition for the mean square of the quark transverse momentum. Of course 
it is different from the naive, gauge dependent definition like {N\sdj_s\N), 
because the physical transverse gluon is participant of this definition. How- 
ever, the expression (p^) is the only possible way to define the {k\)s in the 
gauge theory like QCD. We believe that such definition is the useful gener- 
alization of the transverse momentum conception for the interactive quark 
system. Let us note that the Lorentz transformation in z direction does not 
affect the transverse directions. Thus, the transverse momentum {k\)s cal- 



culated from eq. (|2^) remains unchanged while we passing from the light 
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cone system to the rest frame system where a quark model suppose to be 
formulated. 

Now, let us come back to our definition for {k\)s- In order to cal- 
culate this matrix element, we use the same trick as before: we reduce our 
original problem of the calculating of a nucleon matrix element to the prob- 
lem of a computing of the corresponding vacuum correlation function (pISl) : 

{N\s{t D^fs\N) ~ , . t / dy{0\T{s{t D^)h,uu{0)}\0). (25) 



To estimate the right hand side of the eq.(P^) we introduce an auxiliary 



vacuum correlation function 

t J dy{0\T{sii I D,)s, Mn(0)}|0) = Cg^,, (26) 

where C is constant. From the definition (^) it is clear that the correlator 
we are interested in can be expressed in terms of the constant C: 

i J dy{0\T{s{i D±fs,uu{0)}\0) = -2C. (27) 

At the same time the constant C is given by the correlation function which 
contains Gf,„ and not a covariant derivative Dn: 

C=\ij dy{Q\T{s{i I 5^)s,mm(0)}|0) = 

J dy{0\T{stgG;,^a,,s,uuiO)}\0), (28) 

where we have used the equation of motion and identity^: 

1 ig A" 

D^D^g^^s = '-i^-i^D^D^s - a^^-[D^, D^]s = -m^s + —(Jf^^Gl—s (29) 



Now we can estimate the unknown vacuum correlator (|2^) exactly in the 
same way as we have done before for the correlation function K, see eq. (111 



Collecting all formulae (|J-|2^) together, we arrive to the following final result 
for the mean value of the momentum carried by the strange quark in a 
nucleon: 

_ {N\s{z D^fs\N) ^ {N\-stga,.G';,,^s\N) 
(Ar|ss|Ar) " A{N\ss\N) 

. . . h0.8GeV^f- ^ 0.33GeV^ (30) 

4(ss) d'' V 'S ' V ; 



■^We neglect the term proportional to in the eq. (|2^) . It can be justified by using 
the estimation (|l7|). 
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where dP denotes the dimension of the operator O. For numerical esti- 
mation we use the standard magnitude for the mixed vacuum condensate 
{siga^yG^^j^^s) = 0.8GeV'^{ss). The obtained numerical value (pOD for {k\)s 
looks very reasonable from the phenomenological point of view. 

We close this section with a few remarks. First, the nonvalence nucleon 
matrix elements can be expressed in terms of vacuum condensates in a very 
nice way. All numerical results obtained in such a way look very reasonable. 
As the second remark, we emphasize that a study of nonvalence nucleon 
matrix elements and an analysis of the QCD vacuum structure is one and the 
same problem. We would like to note also, that the nucleon matrix element 
(pop might be very important in the analysis of neutron dipole moment. 
This observation is based on the fact that the so called chromoelectric dipole 
moment of the s quark, related to the operator sg'j^a^uG'l^^^s, in many 
models gets a large factor ~ mg/mg ~ 20 in comparison with a similar d 



quark contribution |T^, [y], At the same time, as we can see from 

(pop there is no any suppression due to the presence of the s quark in the 
corresponding nucleon matrix elements. 



3 Strangeness in the nucleon, channel. 
3.1 Singlet axial constant g^^. 

In this section we discuss the contribution of the strange quarks into the 
nucleon matrix elements similar to eq. (P^D, with the only difference that 
we switch the scalar channel ss into the pseudoscalar one sij^s. In our pre- 
vious study of the scalar channel we concluded that the considerable part 
of the nucleon mass (about 40%) is due to the strange quarkQ. We made 
this estimation by using two following facts: First, we knew the mass of the 
nucleon (left hand side of the eq. (p2D), which is considered as an experi- 
mental data. The second, we calculated independently the matrix element 
(A^|ss|A^). Comparing this theoretical result (p!8| ) with (p2D, we have made 
aforementioned conclusion about a serious deviation from the chiral SU{3) 
limit . 

We want to repeat all these steps for the pseudoscalar channel also. In 
this case the equation analogous to (^) looks as follows: 

2m^°p^75P = +{N\J22mgqt^,q\N) + {N\^G,M~,,\N), (31) 

q 4 TT 

where sum over q is the sum over all light qurks u, d, s and g'^ is the nucleon 
axial constant in the flavor singlet channel. The world average is: g\ = 
0.27 ±0.04, §. 

^In the chiral hmit — > 0, the correspondmg contribution is zero, of course. 
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Now, we would like to repeat all steps which would bring us to the con- 
clusion similar to eq.(^) for the pseudoscalar case. We shall try to answer on 
the following question: what is the strange quark contribution in the formula 
(Pip? Let us recall that in the chiral limit m„ = ma = mg = the nonzero 
contribution comes exclusively from the gluon term in the close analogy with 
formula (0): 

3 

2mg\pi'^^p = -{N\—G^uG^u\N), niu = rud = ms = 0. (32) 
4 vr 

Thus, in order to answer on the question formulated above, we have to esti- 
mate the matrix element 



{p\2mssij5s\p) (33) 

in somewhat independent way Q. First of all, the relevant contribution with 
octet quantum numbers (77) can be easily evaluated by the standard technics. 
One should take the derivative from the octet, anomaly- free, current ~ 
M7m75^ + dltil^d - 2s7^75S. The result is: 

{p\2mssi'^t>s\p)y^ = -m{3F - D)pi-f5p, (34) 

where D ^ O.GSqa and F ~ 0.37gA are the standard SU{3) parameters. One 
could expect that the similar contribution with singlet quantum numbers {r]') 
is also large, although it is zero in the chiral limit where rris = 0. 

We shall estimate the corresponding contribution with 77'- quantum num- 
bers by using our previous trick (|16|). Namely, we reduce our original problem 
of calculation of the nucleon matrix element to the problem of the computa- 
tion of certain vacuum correlation function where we should limit ourself by 
calculating the contribution with singlet quantum numbers only: 

\p) ^ -^^^ J dy{0\T{2m,-sz^,s,m^M0)m. (35) 

In order to make the corresponding estimations we need to know the fol- 
lowing 7]'- matrix elements: (0|sz75s|?7') and {0\ui'j5u\rj') . The PCAC does 
not provide us with the corresponding information, however a quark model 
prejudice suggests that 

{0\-^{ui-f5U - di-f5d)\'JT) ~ {0\-^{uij5U + di-f^d - 2si-f5s)\r]) ~ 



1 - 2(gg) 

{0\ui-f5s\K) ~ (0|^(m75M + di-f^^d + si-f5s)\r]') ~ — (36) 

V 3 Jw 



'We neglect m, d qurk contributions into the formula (pll) by the obvious reasons 
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The strong support in favor that the relations (pq) are to be correct, comes 
from the analysis of the two photon decays of vr, rj, rj', see e.g. All of these 
decay amplitudes have the same Lorentz structure and determined by the 
matrix elements (PBD, therefore, the quark model prediction is found to work 
surprisingly well in this particular case. Combine the formulae (^4|-^) we 
arrive to the following estimation: 

{p\2mssi-f5s\p) = 2mpii75p{-i(3F - D) - } 

^ (—0.3 + 0.16)2mpi75p ~ (— 0.14)2mpi75p, (37) 

where we used 3F — D ~ 0.6 for the numerical estimation. The two terms 
in this formula are the octet and singlet contributions correspondingly. One 
should note, that in spite of the fact that the singlet term is parametrically 
suppressed in the limit rris = 0, this contribution numerically is not small. It 
is only by a factor of two less than the parametrically leading term. 

Now, let us come back to eq.(^). We would like to answer on the previ- 
ously formulated question: what is the s quark contribution into the formula 
(^)? From our estimation (|37| ) we suggest the following pattern of saturation 



of the experimental data for = 0.27 ± 0.04: 

2mp275p(0.27 ± 0.04) = +{p\2m,sij5s\p) + -{p\—G^,G^,\p) 

4 TT 

^ (— 0.14)2mp275p + (+0.41)2mpz75p. (3^ 



Here the first term is due to the strange quark contribution from (pTf ) and the 
second one is due to the gluon contribution. We assign an average number 
0.41 for the gluon contribution in order to match the experimental data for 



9% 



Two remarks are in order. First, The strange-quark and the gluon terms 
contribute with the opposite signs into g^. In the formula for mass (|22|) the 



similar terms interfere constructively, with the same signs. It is very easy 
to understand the difference: in the pseudoscalar channel we have the Gold- 



stone boson, T], whose total contribution is zero into the sum (|3T|) because 
of the octet origin of the rj meson. However, the rj meson contributions into 
the matrix elements {p\2mssi'yr^s\p) and into the gluon operator, taken sepa- 
rately, are not zero. Moreover, its contribution to the {p\2mssi'y5s\p) has the 
opposite sign to the rj' contribution (because of the difference in the quark 
context, see (^)). Even more, it has a parametrical enhancement. We have 
nothing like that in the scalar channel (P^D, where the flavor -singlet states 
dominate. 

The second remark is the observation that, like in the scalar channel, the 
strange quark operator gives a noticeable contribution into the final formula 
(^) in spite of the fact that in the chiral limit the corresponding contribution 
is zero as we mentioned earlier 
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3.2 Singlet axial constant of heavy particles. 

In this section we try to answer on the following question: How one can check 
the estimation ( |38D about noticeable contribution of the strange quark? More 
specific question we would like to know: Is it possible to measure a nucleon 
matrix element where some independent combination of those operators en- 
ters? If answer were "yes" , we would be able to find the contribution of each 
term separately. 

To answer on this question we suggest to consider the weak neutral current 
containing an isoscalar axial component associated with nonvalence quarks 
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{p\cl^,lbC - s-i^-i^s + t7^75t - h-i^-i^h\p) = g^r''^p-i^,-i-,p (39) 

Before to go into details, let us mention, that on the quantum level the 
current divergence of the massive quark field has the following form: 

— — OL ~ 

<9;,Q7/^75Q = 2mQQi-f5Q + -^G^^G^y, (40) 

where the first term is the standard one and the second term is due to the 
anomaly. There are many ways to understand the origin of the anomaly; 
basically it arises from the necessity of the ultraviolet regularization of the 
theory. In the heavy quark mass limit, one can expand 2mQQi'-j^Q term in 
the eq.dO) with the following result ||,||]: 

as ~ GGG 1 
•-G,.G,. + c^ + 0(^ 



2mQQi-i^Q\mQ~.oo = - j^G^uGf^u + c — — + 0(— ) + ... (41) 



where all coefficients, in principle, can be calculated. We are interested, 
however, in the leading term ~ G^yG^j^u only. 



One can easily note that the leading term in the expansion (^l]) has the 
same structure as an anomaly term (|40|) and it goes with the opposite signQ. 
Thus, the terms, which do not depend on mass are canceled out and we left 
with a term which vanishes in the limit niQ — > oo. Such a vanishing of 

Q 

the heavy quark contribution into the nucleon matrix element is in a perfect 
agreement with a physical intuition that a nucleon does not contain any 
heavy quark fields, at least in the limit mq —>■ oo. 

The situation with strange s quark in the formula (^) is much more 
complicated. This quark is not heavy enough to apply the arguments given 

^ The opposite signs of those contributions can be easily understood in terms of the 
PauH-Villars regulator fields with mass Mpy —^oo. As is known these fields are introduced 
into the theory for the regularization purposes and they play crucial role in the calculation 
of the anomaly (^). Regulator contribution is obtained, by definition, by a replacement 
'^Q ^ MpY in the corresponding formula. It goes, by definition, with relative sign minus. 
From such a calculation it is clear that the leading terms which do not depend on mass, 
are canceled out in a full agreement with an explicit formula (^,^l]). 
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above. Thus, we should keep all operators in the formula (1^) for the current 
divergence in the original form: 

{p\2msSi-f5S + ^G/.i/G^z.b) = -g'^''''^2mpi-f5p, (42) 

where we have neglected the terms ~ for c, b, t quarks in according with 

our previous discussion^. 

As we already mentioned, the measurement of the constant g^""^^ (0) 
gives an independent information complementary to the singlet axial constant 
measurement g% (0). If we knew those constants with high enough preci- 
sion, we would be able to find out both nucleon matrix elements: (p|2m5si75s|p) 
and {p\^G ^j,G fj,y\p) . At the moment the experimental errors for g'^"'^^ = 
0.15 ± 0.09, are large: the result is only two standard deviations from 
zero. 

The best we can do at the moment is to estimate g^""^^ from our previous 
calculations (pH]). If we literally take the values —0.14 and 0.41 from the 
formula (^), we get the result for g^""^^ which is compatible with zero: 

- g^^''''^2mpi'-f^p = {p\2msSi-f5S + ^Gf,^Gf,,y\p) ~ 



~ (-0.14 + -0.41)2mp^75p, c/^'™^ ^ 0. (43) 
3 

From our point of view this is an interesting observation which essentially 
says that the strange quark operator together with its anomalous part 
gives nearly vanishing contribution into the nucleon matrix element. As we 
discussed earlier, this is certainly true for any heavy quark. What is surprised 
us, that estimation (^) apperently says that this is true even for s quark 
(which is by no means can be considered as a heavy quark). If we accept this 
point, we should interpret a nonzero magnitude of g'^ (^8|) as a contribution 
coming exclusively from the light u, d quarks and their anomalous parts. 
As we mentioned, the s quark term together with its anomalous part gives 
almost vanishing contribution into eq. (pSj) . Such an interpretation is in a 
very good agreement with the valence quark model philosophy, where the s 
quark does not play any essential role. 

Let us note that this interpretation is very different from the old simplest 
assumption on the spin of the strange quark in the nucleon, see e.g.|^,[^]. In 
our interpretation we understand the strange quark contribution as a joined 
contribution of s field as well as its regulator field (or what is the same, its 
anomalous contribution). 

^ The corresponding estimations even for the hghtest heavy c-quark support this view- 
point, see the next section. 
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4 The charmed quark in the nucleon. 



In this section we would like to extend our analysis for the c quark. The 
reason to do so is twofold: First, the c quark is a heavy enough to use the 
standard 1/mc expansion similar to (|4ll) . Secondly, the charmed quark is 
light enough to get a reasonably large effect from this expansion. 

We start from the pseudoscalar channel and keep only the first term in 
the heavy quark expansion (^Tj): 

Qs ~ —0.41 

(p|d75c|p) ~ -{p\—^Gf,^Gf,^\p) 2mpi'^^p Q.lpi'^^p (44) 

Smcvr QrUc 

where, for the numerical estimate, we use the value (^) for the gluon matrix 
element over nucleon, and nic — l.SGeV for the charmed quark massQ This 
value should be compared with the similar matrix element of the strange 
quark over nucleon: 

-0.14 

(p|si75s|p) ~ — 2mpi'j^p ~ —O.Spi'j^p, (45) 

where we use formula ( |H7| ) for the numerical estimation of the matrix element 
(p|si75s|p). The ratio of these values are in remarkable agreement with the 
ratio of their mass: ^ ~ ~ 7.5. 

ms 0.175GeV 

Our next example is the scalar matrix element. In this case one can 



use the heavy quark expansion similar to formula (|4l|) , but for the scalar 
channel 0: 

{p\cc\p) ^ -ipl^^GlOlJp) ~ ^ • 520MeV^ ^ 



For the numerical estimation in this formula we adopted the value (p3D for 
the gluon matrix element over nucleon. The magnitude (^) for the charmed 
quark is approximately hundred times less than the corresponding matrix 
element for the strange quark ([T8|) : 

{p\ss\p) - 2-4. (47) 

This is in a big contrast with pseudoscalar channel, where the corresponding 
ratio was about a factor ten larger. 

We conclude this section with few remarks. First of all, the matrix ele- 



ments (@),(0) for the charmed quark are expressed in terms of the gluon 
operators. For the heavy quark this is exact consequence of the QCD. Cor- 
rections to these formulae can be easily estimated. One can show that they 



^ Let us note, that this value for mass corresponds to the high enough normaUzation 
point of order of rric- In principle, one should renormalize this value to the low normaliza- 
tion point. We neglect this small logarithmic effect in this paper. 
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are small for the c quark. The problem of the evaluation of the gluon ma- 
trix elements over nucleon is the different problem. However, we believe 
that from the measurements of g'^ (^) and from vrA^ scattering we know 
those matrix elements with a reasonable accuracy. Thus, we expect the same 
accuracy for the matrix elements {p\cc\p) and (p|a75c|p). 

Our next remark is the observation that results (0), (|46|) essentially give 
a normalization for the intrinsic charm quark component in the proton. This 
is very important characteristic of the nucleon. It might play an essential 
role in the explanation of a discrepancy between charm hadroproduction and 
perturbative QCD calculations. We refer to the original paper f\ on 



this subject where the hypothesis of intrinsic charm quarks in the proton 
was introduced. The experimental fit in the framework of this paper [p^ 



suggests that the probability to have an intrinsic charm in the proton is about 
~ 0.3%[^] and (0.86±0.60)%[Q. These numbers can not be related directly 
to the matrix element (^6]) we calculated. However, they give some general 
scale of this phenomenon. We hope that in future, some more sophisticated, 
QCD-based methods, will lead us to deeper understanding of the effects 
related to intrinsic charm component in the nucleon. 



5 Conclusion 

We believe that the main result of the present analysis is the observation 
that non-valence quarks play an important role in the physics of nucleon. 
However we should stress that such an interpretation does not contradict to 
the bag model where the nucleon matrix element 

(Ar|ss|iV) ~ -{ss)V (48) 

is related to some vacuum characteristics (like condensate or volume of the 
bag V) of the model. It is clear that the chiral vacuum condensate of the 
strange quark is large. So, there is no reason to expect that the corresponding 
nucleon matrix element is small. The same argument can be applied for the 
arbitrary mixed vacuum condensates also (they are presumably not zero||17||). 
This information can be translated, in according to (|30|), into the knowledge 
about transverse momentum distribution of the strange quark in a nucleon. 

In our approach this relation vacuum <^=^ nucleon is clearly seen. Thus, 
by studying the vacuum properties of the QCD, we essentially study some 
interesting nucleon matrix elements which can be experimentally measured. 
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